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We study the electronic states of quasi-one-dimensional organic conductors using the single 
band Hubbard model at half-filling. We treat the effects of the on-site Coulomb interaction 
by the fluctuation-exchange (FLEX) method, and calculate the phase diagram and physical 
properties. The calculated pressure dependence of the Neel temperature coincides well with 
the experimental one. We also show that a pseudogap is formed in the density of states near 
the chemical potential and that d-wave superconductivity appears next to the antiferromagnetic 
state. Moreover the NMR relaxation rate increases on cooling in the low-temperature region 
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It is known that the quasi-one-dimcnsional organic 

' superconductors (TMTSF)2X and similar compounds 
have very rich varieties of phases.Ba-' For example, 
(TMTTF)2PF6 is in the spin-Peierls phase at ambi- 
ent pressure. The antiferromagnetic insulating phase 
is stabilized with applied pressure larger than about 
10 kbar. The antiferromagnetic transition tempera- 

\ ture, Tn, becomes maximum 20 K) around 15 kbar. 
(TMTSF)2PF6 is in the antiferromagnetic insulating 

\ phase at ambient pressure. The antiferromagnetism dis- 
appears with increasing pressure to give rise to super- 

' conductivity. The maximum superconducting transition 
temperature is about 1 K. The superconducting phase is 
next to the antiferromagnetic one and it indicates that 
the spin fluctuation is closely connected with the emer- 
gence of the superconductivity. It is also very interest- 
ing that the pseudogap was observed in the density of 
states (lOpS) near the chemical potential in XPS exper- 
iments. lIq' 

[ In the previous theoretical studies, the Hubbard model 
. was mainly studied as the simplest model in these sys- 
' tems. The SDW transition temperature and phase 
i boundary between the SDW and the superconducting 
' phases wemjdiscussed within the mean field, or the RPA 
treatment Shimahara discussed the.spin fluctuations 
as an origin of the superconductivityE-' This work and 
the QMC study indicated that the superconductivity was 
realized and the superconducting order parameter had 
line-nodes on the Fermi surface .EP Dimensional crossover 
was_discussed based on the renormalization group the- 
ory.H^ 

The aim of this letter is to discuss the origin of pseu- 
dogap behavior detected in the XPS experiments, and 
the origin of superconductivity on the basis of the Fermi 
liquid theory, and to understand the phase diagrams of 
TMTSF and TMTTF salts. We show that these phys- 
ical behaviors can be explained naturally in terms of 
spin-fluctuation theory. This is the first study on the 



quasi-onc-dimensional organic conductors where the ef- 
fects of spin fluctuations are taken into account in a self- 
consistent way and where the antiferromagnetism and 
the superconductivity are treated on the same footing. 

These systems have a half hole per TMTSF/TMTTF 
molecule in the HOMO level—and the one-dimensional 
chain is connected weakly.EJ'El^ The dimerization, 
though it is weak, exists and the antibondingdjand which 
crosses the chemical potential is half-fillcd.E3' Therefore 
we consider a single-band Hubbard model whose band 
dispersion is given as 

e{kx,ky) — 2to cos ak^ + 2t i cos bky -f 2t2 cos 2bky , ( 1 ) 

where |to| ^ \ti\ ^ 1^2 1 and a and b arc the effective lat- 
tice constants. It was known that the cos 2bky term was 
necessary to appropriately reproduce the shape of the 
Fermi surface.© The physical role of this term is that it 
breaks the nesting condition, in other words, that it in- 
troduces frustration. It was also known that t2 was pro- 
portional to ti/to, when this term was derived from the 
quarter-filled Hubbard model. We study the single-band 
Hubbard model at half-filling and assume the screening 
of long range Coulomb interactions for simplicity. This is 
justified because the antiferromagnetic fluctuations play 
the essential roles. 

In this study, we put to — —1, t2 — 0.8i^/io and the 
on-site Coulomb interaction U = 3.0, and vary ti{< 0) 
as a parameter. Here \t2 \ is rather large compared to tb£ 
model derived from the quarter-filled Hubbard modelE' 
Such a choice is mainly due to the convenience of the 
numerical calculation because we use a finite k-mesh. 
We implicitly use 128x128 k-mesh and 256 Matsubara 
frequencies in this study. 

To tackle the present model we employ the fluctua- 
tion exchange (FLEX) method, which is a kind of self- 
consistent perturbation theory with respect to U. The 
FLEX method has been considered to be advantageous 
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for systems with large spin fluctuations. It was known 
that this method gave the Green function which agreed 
with the QMC oneJor the square-lattice Hubbard model 
with moderate US3 This method was applied to the 
study of high-Tc cuprates, superconducting latMeLcom- 
pounds and k-(BEDT-TTF) salts and so onJlfl'BEJ and 
predicted the possible d-wave supercondiictivity. The 
QMC study also supported this result.ErliZP 

Here, we explain the FLEX method. The Dyson equa- 
tion is written as 

{G(k,e„)}-i = {GO(k, -S(k,e„), (2) 

where G'^(k, e) is the unperturbed Green function. The 
self-energy at temperature T is given by 

S(k, e„) = G(k - q, e„ - iui) ■ 

q,; 

X Qx(-Hq,c.o + ^x<+Hq,'^0-x°(q,^o) , (3) 

X(±)(q,c.0 = x°(q,c.O • {1 ± C/x°(q,c.O}"' , (4) 
X°(q,c^O = -r^G(q + k,c^; + e„)G(k,e„), (5) 

k,n 

where e„ — {2n + l)7rT andtj; = 2lTrT. Equations (2)-(5) 
are calculated together with the equation for the chem- 
ical potential ^ given by = 2T „ G(p, e„)e""° , 
where 1 in this system. 

To determine the magnetic transition temperature Tn, 
we calculate the Stoner factor without vertex corrections, 
as — maxjj. I U ■ x^0^i^ = 0) }• The antiferromagnetic 
critical points are determined by the Stoner criterion, 
as = 1- 

We also determine the superconducting transition 
temperature Tc by solving the linearized Eliashberg 
equation with respect to the singlet-pairing order param- 
eter, 0(-k, e„) = -|-(/)(k, e„), 

A • 0(k, e„) ^-TJ2 - q, e« - em) 
q,m 

xG(q, e„)G(-q, -e„) • (/)(q, e„), (6) 
V(k,uji) = ^U\^-Hk,L,i) - ^U\^+\-k,LUi) + U, (7) 

where Tc is given by the condition A = 1. 

The theories by Mermin and Wagner, and Hohen- 
berg prohibit finite Tn and Tc in two dimensions. It 
was well known that as satisfies this condition, because 
FLEX treats the spin fluctuations self-consistently.H 
We estimate Tn by the condition as = aN, where 
(1 — aN)^^ ~ O(IOO). This implies that the weak mag- 
netic coupling between layers J± makes the system or- 
dered three-dimensionally. On the other hand, A = 1 is 
fulfilled at finite Tc using eq. (6). However the Tc given 
by the Eliashberg equation is reliable in many cases. 

First, let us examine the phase diagram on the plane of 
|ti| and T. Each line corresponds to Tn for different aN. 
With increasing \ti\ { < 0.3) Tn increases gradually and 
has a broad peak (Tn ~ 0.4) around |ti| = 0.3. On the 
other hand, Tn decreases rather rapidly with increasing 
1^1 1 ( ^ 0.4). This behavior is consistent with experi- 



ments if one regards \ti \ as applied pressure. We do not 
have a plot for \ti\ < 0.1, where k-mesh and Matsubara 
frequencies are insufficient to determine Tn numerically. 

An inset of Fig. 1 shows the phase diagram at U = 
1.8 where x(q, tj) is calculated using G°(k, lj) and the 
Neel temperature is estimated by the Stoner criterion, 
aN = 1- The value of U is chosen in order to make Tn 
drop at |ii| ^ 0.4. In this case Tn does not descend with 
decreasing \ti\. Therefore the decrease of Tn for small 
in Fig. 1 is due to the self-consistency. We note that 
the estimated Tn using G"(k, a;) is about 10 times larger 
than that using self-consistent G(k, lj), at ti = —0.3 and 
U = 3.0. 

Fig. 1. The Neel temperature, Tjm, determined self-consistently 
as a function of \ti\ at U = 3.0, and on = 0.992 (square) and 
On = 0.994 (circle). An inset shows the Tn using in x a-t 
U = 1.8. The incommensurate region corresponds to the closed 
symbols. 



In this model, the nesting vector Q is commensurate, 
or (7r,7r), except in the large |ti| and small T region. 
In Fig. 1 the incommensurate region corresponds to the 
closed circles and squares for T ;< 0.02. The nesting vec- 
tor is incommensurate in the ky direction, but on the 
other hand it is always commensurate in the direction 
whereas the RPA result is notfl't^ In other words, the 
SDW phase is locked to n in the kx direction. This pSjiljt 
is consistent with the experimental nesting vectoil23'ElP 
and is first realized in the self-consistent calculation. We 
note that the nesting vector depends on the details of the 
Fermi surface. Therefore it is necessary to use a realistic 
Fermi surface in order to compare the calculated nesting 
vector with the experimental one for all regions of the 
phase diagram. 

Next let us see the Fermi surface. The Fermi sur- 
face in the present model is flatter near kx ^ ±7r and 
pointed near kx ^ due to the cos{2bky) term. In Fig. 2, 
the Fermi surface is shown at ti = —0.4, and T = 0.05 
and U = 3.0. The Fermi surface with interaction has a 
smaller band dispersion in the ky direction, at the same 
time, the nesting condition is better than the one with- 
out interaction. 

Fig. 2. Fermi surfaces at U = 3.0 (a solid line) and U = {a 
broken line) at ti = -0.4 and T = 0.05. 



We plot the DOS in Fig. 3. The pseudogap emerges 
near the chemical potential due to the spin fluctuations 
and it evolves with decreasing |ti | and T. This result has 
close connection with the pseudogap near the chemif;a.l 
potential found in the XPS and ARPES experiments BeI' 
We note that this model is effective near the chemical 
potential. 

We also show the uniform spin susceptibility 
(X" (0, uj = 0)) and the NMR relaxation rate ((TiT)"! = 
lim^j^o Im {k . iS) / ijj) in Fig. 4. With decreasing 
temperature, the uniform spin susceptibility decreases 
gradually and this behavior reflects the evolution of 
the pseudogap in the DOS near the chemical poten- 
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Fig. 3. (a)The density of states at ti = —0.4 and U = 3.0, and 
T = 0.05 (a solid line) and T = 0.1 (a broken line), and at 
ti = —0.4 and C/ = (a dotted line). (b)Tlie density of states at 
U = 3.0 and T = 0.1 and \ti\ =0.1, 0.2, 0.3, and 0.4, respectively. 



tial. (TiT)^^ increases with decreasing temperature 
in the low-temperature region. Staggered susceptibility 
(X (Qj'^ = 0)) is also the case. These results seem to 
be consistent with experimental results. LaL^ 

Fig. 4. The temperature dependence of x'^Ho,'^ = 0) and 
1/TiT at ii = -0.4 and U = 3.0. 



One can calculate the Green function at lower temper- 
ature than T ~ 0.01, though the 256 Matsubara frequen- 
cies are not sufficient. The superconducting condition is 
fulfilled at ti = -0.45 and T < 0.007 in this model, when 
we use 128x128 k-mesh and 512 Matsubara frequencies. 
We plot (j){kF,0) corresponding to the singlet SC order 
parameter in Fig. 5(a), where T = 0.007 and A in eq. (6) 
equals 1.02. One can clearly see that nodes exist and a 
c?-wave-like superconductivity is realized. This is consis- 
tent with the NMR experiment, jSjzliich predicts the SC 
order parameter with line-nodesO-* There exists a dip 
at fcy ~ and a peak at ky ^ tt. Figure 5(b) shows 
a contour plot of 0(k, — 0). At ky — 0, the peak in 
(/>(k, w = 0) is closer to k^ — 7r/2 than the Fermi surface, 
and at ky — ±7r the situation is reversed. These dip and 
peak structures depend on the manner of warping of the 
Fermi surface and disappear when 1^2 1 is small. We have 
checked that A for the triplet case is about one third 
when the singlet SC is fulfilled. Therefore the triplet 
superconductivity is not realized in this model. 

Fig. 5. (a)Tlie superconducting order parameter on the Fermi 
surface, (f>{k-p,i^ = 0), from ky = — tt to ky = tt and (b)the Fermi 
surface and the contour plot of (f>(\s.,u} = 0), at ti = —0.45, 
U = 3.0 and T = 0.07. 



Below we focus on the obtained phase diagram in 
Fig.l. The self-consistent Tn declines with decreasing 
|ti|, whereas the non-self-consistent Tn does not descend 
with decreasing \ti\. At the same time, the DOS at the 
chemical potential decreases with decreasing |ti| because 
of the strong spin fluctuations in the self-consistent cal- 
culation as shown in Fig. 3. This effect is not taken 
into account when one uses These two behaviors 
are reasonably understood through the Kramers-Kronig 
relation, 



X 



0) = - P 



duj, (8) 



where the main contribution comes from small uj. 
The integrand can be written as l/uj Imx*^^-'(k, cj) ~ 
l/a;(Imx°(k, Lj))(l — J7x°(k, o;))""^ for large as, and 



l/uj Imx*'(k, +iO) at zero temperature is estimated as 
lim -Imx"(k, +iO) ^ V ImG(q, 0) ImG(q-|-k, 0). (9) 

The right-hand side of eq. (^ is roughly proportional 
to the square of the DOS at the Fermi level. Therefore 
x'^^k, 0) decreases together with the evolution of the 
pseudogap in the DOS near the Fermi level. This is the 
main reason why Tn decreases for small \ti\ in Fig. 1. 

One can interpret the reduction of Tn for small \ti\ 
in Fig. 1 in connection with the 5=1/2 coupled chain 
Heisenberg model.ES' The magnetization at T = de- 
creases due to the zero-point fluctuation in decreasing 
the interchain coupling. 

On the other hand, when \ti \ is large enough, then \t2\ 
is also large, and the nesting condition itself determines 
Tn. Hence Tn drops rapidly as seen in the RPA and 
mean field calculations. 

The eff'ect of applied pressure can be considered to 
increase the transverse transfer integrals and the band 
width. But it is sufficient only to enlarge the transverse 
transfer integrals in order to discuss the overall phase 
diagram in TMTTF/TMTSF salts. Experimentally, the 
Neel temperature (or SDW transition temperature) in- 
creases with increasing pressure, whereas it decreases 
and finally drops rapidly near the phase boundary with 
the superconducting phase.© We have derived, for the 
first time, the phase diagram which coincides with the 
experimental one by considering the effects of the strong 
spin fluctuations self-consistently. 

Next we compare the present results quantatively with 
the experimental results. If |io| is set to 1000 K, T = 

0. 01 correspond to 10 K. These results do not correspond 
well with the experimental ones because U and 1^2 1 are 
rather large in the present study due to the numerical 
calculations. But Tn is reduced with decreasing U and 
with increasing |<2|, and the value of ti at which Tn takes 
its maximum also behaves in a similar way. Then one can 
adjust U and \t2\ such that Tn and ti are close to the 
experimental ones. 

We have calculated a quasi-one-dimensional model 
whose lattice has zigzag hopping parameters, and have 
shown that the results are almost the same as that of the 
current model. We have also calculated the quarter-filled 
Hubbard model without dimerization. The temperature 
dependence of the Stoner factor as a function of ti is 
very similar to the case of the half-filled one in Fig. 1, 

1. e., Tn has a peak as a function of ti. But the resulting 
maximum Tn seems to be less than about 5 x 10^'^ and 
the dimerization takes an important role quantitatively. 
These results will be published elsewhere. 

The recent experiments revealed that the antiferro- 
magnetic phas£s were accompanied by charge dispro- 
portionationsH Naturally the importajice-of long-range 
Coulomb interactions was pointed out.ca'EZi' The corre- 
sponding degrees of freedom, which are dropped in the 
current model at half-filling, will not be negligible be- 
cause of the weak dimerization. Further studies includ- 
ing this effect are necessary in connection with the nest- 
ing vector. 
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In summary, we have studied the model Hamihonian 
of the quasi-one-dimensional organic conductors by the 
FLEX method. We have successfully explained the pres- 
sure dependence of the Neel temperature as a result of 
the evolution of the pseudogap in the DOS near the 
chemical potential, and the NMR relaxation rate which 
increases on cooling in the low-temperature region. We 
have also found a d-wave-like superconductivity with 
line-nodes on the Fermi surface. 
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